DEPARTMENT OF MATHEMATICS & STATISTICS 
COLLEGE OF ENGINEERING AND TECHNOLOGY 
OLD DOMINION UNIVERSITY 
NORFOLK, VIRGINIA 23529 


THE REDUCTION OF DUCTED FAN ENGINE NOISE VIA A 
BOUNDARY INTEGRAL EQUATION METHOD 

By 

Principal Investigator: 

Dr. J Tweed 

and 

Co-Principal Investigator. 

Dr. M. Dunn 

Final Report 

For the period April 1996 through April 1997 

Prepared for 

National Aeronautics and Space Administration 

Lewis Research Center 

Attn.: Dennis Huff, Technical Officer 

Mail Stop 77-6 

21000 Brookpark Road 

Cleveland, OH 44135 


Under 

Grant Number NCC-3-469 
ODURF #136161 


February 1997 




Final Report: April 1996 - April 1997 

The Reduction of Ducted Fan Engine Noise Via a Boundary Integral Equation Method 


M.H. Dunn and J. Tweed 
Department of Mathematics and Statistics 
Old Dominion University 


Introduction 

The development of a Boundary Integral Equation Method (BEM) for the prediction of 
ducted fan engine noise is discussed. The method is motivated by the need for an efficient and 
versatile computational tool to assist in parametric noise reduction studies. In this research, the 
work in reference 1 was extended to include passive noise control treatment on the duct interior 

The BEM considers the scattering of incident sound generated by spinning point thrust 
dipoles in a uniform flow field by a thin cylindrical duct. The acoustic field is written as a 
superposition of spinning modes. Modal coefficients of acoustic pressure are calculated term by 
term. 

The BEM theoretical framework is based on Helmholtz potential theory. A boundary value 
problem is converted to a boundary integral equation formulation with unknown single and double 
layer densities on the duct wall After solving for the unknown densities, the acoustic field is easily 
calculated. 

The main feature of the BIEM is the ability to compute any portion of the sound field without 
the need to compute the entire field Other noise prediction methods such as CFD and Finite 
Element methods lack this property Additional BEEM attributes include versatility, ease of use, 
rapid noise predictions, coupling of propagation and radiation both forward and aft, 
implementable on midrange personal computers, and valid over a wide range of frequencies 




Accomplishments 


Given a locally reacting, segmented liner model on the duct interior, a system of boundary 
integro-differential equations for the scattered acoustic pressure field was derived (references 2 
and 3). The equations are hypersingular and require special analytical and computational 
techniques to solve. These methods are near completion. 

For small inflow Mach numbers (M<0.4), the complexity of the integra equations can be 
significantly reduced. The small Mach number BEEM has been determined to be valid at takeoff 
and approach. Mathematical details of the derivation and solution techniques will in reference 2 
BEEM methodology and several passive and active noise control studies were presented at 
the NASA AST Engine/Nacelle Noise Workshop (reference 3). In this presentation, the accuracy, 
versatility and simplicity of the BEM were demonstrated. 

A computer program for predicting ducted fan engine noise (TBEM3D) and user manual 
(reference 4) has been written. The TBEM3D code was designed for general use by the 
aeroacoustics community and is available electronically. 

BEEM methods were applied to sound radiation and propagation in two dimensions. The 
resulting methodology is valid for the shielding of sound by thin strips (wings) and the radiation 
and propagation of sound in a finite length channel. The 2-D channel configuration is an excellent 
tool for studying the effects of a scarfed inlet on noise directivity (see reference 5) A general use 
2D noise prediction computer program (TBIEM2D) and user manual are near completion 




Conclusions 


The effectiveness of BIEM as a tool for active and passive noise control has been clearly 
demonstrated. Due to the simplicity, speed, accuracy, and versatility of BEEM, TBEM3D is far 
superior to wave envelope and CFD codes for conducting parametric noise reduction studies. 

The implicit coupling of sound radiation and propagation both forward and aft are unique 
features of BIEM In addition, the far field radiation condition is satisfied implicitly 
Consequently, special boundary conditions at an artificial far field boundary and at the duct s inlet 
and exhaust planes are not required These theoretical niceties further enhance the simplicity of 
BEM. 
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Ab stract 

A computationally efficient Boundary Integral 
Equation Method (BIEM) for the prediction of ducted 
fan engine noise is presented. The key features of the 
BIEM are its versatility and the ability to compute rap- 
idly any portion of the sound field without the need to 
compute the entire field. Governing equations for the 
BIEM are based on the assumptions that all acoustic 
processes are linear, generate spinning modes, and oc- 
cur in a uniform flow field. An exterior boundary value 
problem (BVP) is defined that describes the scattering 
of incident sound by an engine duct with arbitrary pro- 
file. Boundary conditions on the duct walls are derived 
that allow for passive noise control treatment The BVP 
is recast as a system of hypersinguiar boundary integral 
equations for the unknown duct surface quantities. 
BIEM solution methodology is demonstrated for the 
scattering of incident sound by a thin cylindrical duct 
with hard walls. Numerical studies are conducted for 
various engine parameters and continuous portions of 
the total pressure field are computed. Radiation and 
duct propagation results obtained are in agreement with 
the classical results of spinning mode theory for infinite 
ducts. 


Notation 

denotes that a quantity is dimensional 
when appearing over a variable 
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stretched, moving frame 
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of acoustic velocity field 
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d/oN normal derivative operator in stretched, 

moving frame (with respect to outward 
facing normal to stretched duct surface) 
surface acoustic impedance for duct 
exterior (interior) 
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Ducted fan engine noise is dominated by the fan 
component at takeoff and approach. Community expo- 
sure to the high leveis of radiated fan noise at these 
conditions is significant. The reduction of tonal noise 
produced by the rotating components of high bypass 
turbofan engines is therefore of primary concern to the 
aeroacoustician. The design of active and passive noise 
abatement technology can be facilitated by advanced 
analytical tools for predicting the radiated sound from 
engine ducts. To be useful in design studies, prediction 
tools should be fast, versatile, accurate, and tmple- 
mentable on mainstream computer systems. The ability 
to compute only a portion of the sound field without the 
need to calculate the entire field is an important attrib- 
ute in conducting rapid noise predictions. Computa- 
tional approaches such as Finite Element Methods 
(FEM) and Computational Aeroacoustics (CAA) meth- 
ods lack this property. For this reason, farfieid noise 
calculations using FEM or CAA require vast amounts of 
computational time and computer storage. Therefore, 
the use of FEM and CAA for parametric studies in noise 
abatement research is limited. 

In this paper, a Boundary Integral Equation Method 
for the prediction of ducted fan engine noise is pre- 
sented. The method is based on the equations of line- 
arized acoustics with uniform inflow. A scattering ap- 
proach is adopted in which the acoustic pressure field is 
split into known incident and unknown scattered com- 
ponents. The source process is assumed to generate an 
incident pressure field that can be represented by a su- 
perposition of spinning modes. In a frame of reference 
moving with the engine duct and in regions of space not 
occupied by acoustic sources or scattering surfaces, the 
components of acoustic pressure are governed by 
Helmholtz' equation. An exterior boundary value 
problem is obtained by the inclusion of boundary con- 
ditions on the duct surfaces. The most general form of 
the boundary conditions allows for a spatially varying, 
locally reacting liner model on the duct surface. 

By considering special values of the specific 
acoustic impedance in the boundary conditions, the 
classical Dirichlet and Neumann boundary values are 
obtained. The boundary value problem is then solved 
by expressing the scattered pressure field in terms of 
double and single layer Helmholtz potentials with un- 
known densities that are related to surface pressure and 
the normal derivative of surface pressure, respectively. 
Application of the boundary conditions to the layer rep- 
resentation yields a system of one-dimensional, hypers- 
mgular boundary integral equations for the unknown 
layer densities. The source terms for the system are 
related to the known incident pressure and its normal 


derivative. This system of boundary integral equations 
and method of solution comprise the BfEM. 

The system of boundary integral equations is valid 
for engine ducts with arbitrary profile. If, however, the 
duct is approximated by an infinitesimally thin cylindri- 
cal tube, the complexity of the integral equation kernels 
is substantially reduced. 

Analytical results will be presented that separate 
the singular and logarithmic portions of the integral 
equation kernels from the bounded parts. This analysis 
is significant because calculations involving singular 
and logarithmic integrals are available in closed form, 
thus avoiding time consuming, customized numerical 
integration techniques. 

To demonstrate the BIEM, the solution procedure 
for a thin pipe geometry with hard wall boundary con- 
ditions is presented. A collection of spinning point di- 
poles located inside the duct are used to simulate the 
loading component of the fan noise and generate the 
incident pressure field. Several sets of engine operating 
parameters are considered in this study. Various re- 
searchers 1 ’ 5 have employed boundary integral tech- 
niques to solve this problem. Differences in the present 
work relative to the referenced works appear in the 
conclusions section of this paper. 

Boundary Vaiue Problem-Derivation 

In the analysis that follows, all quantities have been 
nondimensionalized; length by r mtt , mass by 
and time by Q" 1 . 

We consider an engine fan surrounded by an axi- 
symmetric, nondeformable duct of arbitrary profile 
translating in the +z (axial) direction with uniform 
speed V F (see figure 1). The fan is composed of N 9 
equally spaced biades and rotates with shaft soeed n. 
The incident acoustic pressure field generated by the fan 
is known. Linear conditions are assumed to apply and 
the inflow is uniform. 



Figure 1: Duct Geometry 
(Cylindrical, Stationary Coordinates) 
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Governing Differential Equations 


The total acoustic pressure in the sound tie id is 
written as a sum of incident and scattered parts. 

P ;(r # M/ > 2 ,t) = p;(r,n/,z ? t) + p;(r l v,z,t) (l) 
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for the m-th coefficient. Using the definitions m (6), 
the momentum equation (3) can be written as 


In regions of space that contain no scattering surfaces, 
p' is governed by the homogeneous wave equation. 
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Total acoustic pressure and acoustic velocity are 
related through the normal component of the acoustic 
momentum equation. 

oujr.v.z.t) 5p[ (r,v|/,z,t) 

+ = 0 (j) 

at an 

In a frame of reference moving with the duct, all 
dependent acoustic variables can be expressed as linear 
superpositions of spinning modes. For example, the 
scattered pressure has the form 


where is the axial component of the outward facing, 
unit normal to the stretched duct surface. Equation (8) 
is valid for M F > 0. • If the duct is stationary (M F = 0 ), 
then (7) and (8) reduce to 


and 




k : 


P,=0 
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1 5P,(r,Z) 
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= 0. 


(9) 

( 10 ) 


„ Equations (7-10) are valid for points not lying on the 

p' (r, \\f ,Z, t) = £ P> m (r,Z)e (4a) surface of the stretched duct. 

m/Nj = 


and the acoustic velocity is written 

u„(r, V ,Z,t)= (4b) 

nvN B *-« 

where the stretched, moving axial coordinate Z is given 

by 


The BIEM calculates modal amplitudes in (4) term by 
term. For notationai convenience, the superscript m on 
the modal coefficients is dropped hereafter. 

Define the dependent variables Q, and <X> N by 

Q J (r,Z) = P 1 (r.Z)e‘ M - z (6a) 

and 

0 N (r,Z) = (J N (r,Z) e‘* M ' z , (6b) 


Boundary Conditions 

To meet noise certification levels it is necessary to 
treat the engine duct with passive noise suppression 
technology. In this work, the duct treatment is modeled 
by a locally reacting, axially varying liner. 

Define the surface functions Qf and by 

Qf (Z) = lim Q,(r,Z) Zs[a,b] (11a) 
r— r sU! 

and 

<t>^(z)= lim <J> N (r,Z) Ze[a,b] (lib) 

r-t-r5(Z) 

with similar definitions for Qf, Qf, normal velocity, 
and all pressure components. 

Myers 6 has shown that if s'(Z) represents the 
specific acoustic surface impedance, then, in the 
stretched, moving frame of reference, the modal coef- 
ficients of velocity and pressure satisfy the boundary 
equation 


with similar definitions for the total and incident pres- 
sures. Combining (4-6) with (2) yields the two dimen- 
sional Helmholtz equation 


/ 

<t>S(Z) + 




mi a)Qr(z) 
P ozjc-U) 


Z e[a,b].(12) 
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Equation (12) represents two boundary conditions, one 
for the extenor surface and the other for the intenor 
surface. The boundary conditions for the stationary 
case reduce to 


Classical boundary conditions are obtained for the 
stationary case. Combining (13) and (19) yields the 
Robin’s boundary conditions 


u=(z) 


p*(z) 

C*(z) 


= 0 


Z s[a, b]. 


(13) 


p,*(z) i ( ap, y = 

£*(z) ick y on J 


0 Z e[a, b]. 


( 20 ) 


By considering special values for the functions 
in (12) and (13), several boundary conditions of interest 
are derivable. Total sound absorption is achieved if 


C,* = C,~ ~ constant — ► 0. (14) 

Total sound reflection (hard walls) arises for 

= C" = constant -> oo. (15) 


Of particular relevance to actual engine ducts, is the 
case of a hard exterior surface (C * — >oo) and lined in- 
terior. The boundary conditions for this situation are 


(z) = 0 Z e[a,b] 


(16a) 


and 




1 + 


Mp d jQr(z) 


3 5Zjc(Z) 


= 0 Zs[a,b], (16b) 


We now evaluate the momentum equation on the 
duct surface to obtain a relation between surface pres- 
sure and surface velocity. Let 


{ 3Q« 

( 5N 


(Z) = 


!im 


5Q,( r.Z) 

JN 


(17) 


Dinchlet or Neumann conditions are obtained by appli- 
cation of (14) or (15), respectively, to (20). 

We complete the derivation of the boundary value 
problem by requiring the dependent acoustic variables 
to satisfy the appropriate Sommerfeld radiation condi- 
tion. Additional conditions may apply depending on the 
smoothness of the duct wails. If there are points on the 
duct that do not possess continuously turning tangent 
planes, then edge conditions specifying the behavior of 
the acoustfc pressure at these points must be provided'. 
Edge behavior is determined from the physics of the 
problem together with an asymptotic analysis of the 
governing equations in a neighborhood of the edge. 
This problem occurs at the leading and trailing edges of 
the thin pipe approximation, for example. 

Summary 

The above analysis describes a uniquely solvable 
two dimensional boundary value problem for the scat- 
tered acoustic pressure in the sound field. For the sta- 
tionary case, the boundary value problem is denned by 
(1), (9), (19), (20), the radiation condition, and the edge 
conditions (if applicable). If the duct is in motion, then 
(1), (6), (7), (12), (18), radiation and edge conditions 
completely define the boundary value problem. 

Boundary Integral Equation Formulation 


with normal derivatives of other dependent acoustic 
variables evaluated on the duct surfaces similarly de- 
fined. Then from (8), (10), and (1 1) we get 



Z c[a, b] V1 F >0 


and 



Zs[a,b] M f = 0 


In this section, we reformulate the boundary value 
problem by deriving one dimensional boundary integral 
equations in which the boundary functions Qf and Qf, 
are unknown, where 

«< z K§r)“< z >- <-» 

Once the scattered boundary functions are determined, 
the scattered pressure in the sound field is calculated 
pointwise via a Helmholtz layer representation that sat- 
isfies the Sommerfeld radiation condition implicitly. 

Helmholtz Laver Representation 

The Green’s function for the two dimensional 
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Helmholtz operator in (7) can be written 


1 f e 

Q{v,r' ,Z-Z') = — cosmvi/-^— d\|/ (22) 


D s [f](z) = J f(Z')D* (Z,Z')dZ' 

a 

r — r5(Z) 


where 


R = ^r z +r ,2 -2rr' cos v|/ + (z - Z')“ . (23) s * [f](Z) = J f (Z')S^ (Z,Z’)dZ* 

a 

s;(z.z’).^j-(z’),.iim. 1 


Using resuits mm Helmholtz potential theory* and 
(21-23), the solution of (7) is expressed as the sum of 
single and double layer Helmholtz potentials with den- 
sities Q* and Q ; , respectively. Thus, 


oN 


Q,(r,Z) = s*[Q;](r,Z)-d"[Q;](r,Z)- 
s-[Q;](r,Z)-d-[Qr](r,Z) 


(24) 


r — ( Z) 


D-;[f](z)=Jf(Z')D^(Z,Z’)dZ' 

a 

°N ( Z, Z') = — -I ■ ( Z') ^lim 2 } , 


( 28 ) 


(29) 


(30) 


Equation (24) is valid everywhere except for(r,Z) on 
the stretched duct surface. The field operators s~ and 
d* are defined by 
s ± [f](r,2) = 

(25) 


1 

2n 


f(Z') 


lim G 

-'o(Z’) 


J ± (Z')d Z' 


and 


d ± [f](r,Z) = 


in 


f(Z') 


r 50 
lim _ , ; 

—r5(z-) SN 


J t (Z')dZ' 


(26) 


where f is some sufficiently smooth function and 
J~[Z')dZ' are the elements of arclength along the 
curves r = r^ ( Z'). 

If the boundary functions were known, then the 
scattered acoustic field could be obtained from (24). 

Singular Boundan/ Ope rator Notation 


S|[f](z) = Jf(Z')Si(Z,Z')dZ' 


s t(z,z') = 4 r i-(z')to X 

r-»r^(Z) 

D|[f](z) = If(Z')Di(Z,Z')dZ' 

a 

D z(Z,Z') = 2^ J ’( Z ') f .i^ z .)az5N' 

r-'£(Z) 


(3D 


(32) 


Note that operators are denoted by bold face type. The 
kernels are singular for Z — Z' = 0. The nature ot the 
singularities is examined below. 

Singular Kern el Analysis 

We list here, without proof, the asymptotic proper- 
ties for |Z-Z'| « 1 of the above kernels. The singular 
character of the kernels is obtained by local analyses, 
the details of which will appear in a future publication. 


In order to apply the boundary conditions to the 
field equation (24), it is necessary to evaluate the 2-D 
layer potentials and their derivatives on the stretched 
duct surface. The resulting i-D boundary operators are 
singular. For Z e[a,b], define the following singular 
boundary operator - kernel pairs: 


S t (Z,Z , )=gfi( z ) ln ! Z_Z 'h K '(Z>Z') (33) 

(34) 




g* (Z)ln|Z-Z'|+K.,(Z,Z') 


s I [f](z) = lf(z')sMz. z ')d z ' 


S r (Z,Z') 


— J r (Z') lim G 

2tT r-'i(Z') 

'-+r$U) 


(27) 
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ILm z ,d 1 [f]( r ,Z) = -l f ( z ) + Di [ f ]( Z ) (40) 


s * (Z.Z') = 


’,h(z) 


Z-Z' 

gj 3 (z)ln|Z-Z'!-K J (Z,Z') 


Dj(Z,Z* 


g;,(z) ] g| t (z) i 

(Z-Z') 2 z-z ' 
gi(Z)ln|Z-Z'| + K.(Z f Z') 


(36) 


ci&)^r si t f ^ r ’ z)= i f(z)+ss t f I (zl (41) 

,iS.sr d ‘l f ] ( '- z )- D 5[ f Kz) (42) 

Boundary Integral Equations 


S;(Z,Z') = 


;z) 


Z-Z’ 

gj 3 (Z) ln|Z -Z')+ K, (Z,Z') 


D;(Z,Z’) = 


: (Z) 


g|,(Z) , ! 

(Z-Z') : z ~ z ' 

g* 3 (z)ln|Z-Z’| + K 6 (Z,Z') 


(37) 


(38) 


where and K- are known continuous functions. 

The leading behavior for the single layer kernel 
(33) is logarithmic. Therefore, the associated operator 
is weakly singular. The leading terms for the kernels 
(34,35,37) are of the Cauchy type. Whilst, the kernels 
(36,38) are of the strongly singular Hadamard type. 
Consequently, integrals involved with the kernels (34- 
38) are divergent and must be interpreted in the finite 
part sense . 

All of the above kernels have the logarithmic por- 
tions extracted. Integrations involving these terms are 
defined but difficult to achieve numerically. This 
problem is mitigated by the development of analytical 
results for the associated operators. Examples of this, 
as well as analytical results for the Cauchy and 
Hadamard terms, are presented in the results section. 

Calculations involving the continuous portions of 
the kernels are performed by straightforward numerical 
integration. 

Jump Relations from Potential Theory 

Using the above operator notation, we state conti- 
nuity properties for the single and double layers as a 
field point approaches the surface from the exterior of 
the duct. For sufficiently smooth f and (except pos- 
sibly at Z = a,b), we have the following results* for 
Z € ( a , b ) . 


A system of integral equations for Q* and Q ^ is 
derived by applying the boundary conditions (12) and 
(18) to (24) and (39-42). 

We begin by deriving some preliminary resuits. 
From (1), (6a), (1 la), (24), and (39-40) we write 


Qf (z) = Qf (z)+s* j(z) 

+(l-D*)[Q?](z) Z e[a,b] 
and 

^(z).^g.( 2 ).s;[Q 5 ](z) 

-^( z )- D z[ Q f]( z ) Z e[a,b] 
Combining (1), (6a), (17), (24), and (41-42) yields 


£ 0 , 

oN 


(Z) = 


fQi. 

<5N 


(Z)- 


(45) 


(l + S^ )[Q£ ](z) + [Q f ]( Z) Z e[a,l 


Define the unknown surface vector functions q * by 

q- (z) = (q* (Z),Q= (Z)) T . (46) 


Combining the boundary conditions (12) and (18) with 
results (33-46) yields the system of integral equations 

e"^ Z 'Kf[<r](Z’)dZ' 

(47) 

+ K?[q*](z) = qf (z) Zs(a,b) 

for the four unknown surface functions. The vector 
function q ( s is known from the incident pressure field 
and the integral operators Kf and Kf have the general 
form 


J 


!im ] S i [f]( r ,Z) = S £ [f](Z) 


(39) 
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= 0 Z e[a, b]. (50) 


A( Z)q* ( Z) + B(Z 


5q*(Z) 

dZ 


c(z) 


f q s (Z') 


z-Z'Y 


-dZ'-D(z) 


q ~ ( z ' 


z-z 


— dZ' +,(48) 


E(Z)Jq s (Z')lniZ-Z'|dZ' + K B [q £ ](Z) 


/A 


gQ, 

<3N 


= Qn± 


fQi 

dr 


Since the incident pressure and its derivatives are con- 
tinuous across the duct surface we add the extenor and 
intenor equations in (50) to get 

Qn + Qn=° Ze[a,b]. (51) 


where A,...,E are matrices of known functions and K B 
is an integral operator with continuous kernels. The 
matrix functions are determined from the coefficients in 
(33-38) and depend on the surface impedances and the 
duct curves r = r^(Z')- Explicit expressions are 
lengthy and will not be presented here. 

Examination of (48) indicates that (47) is a system 
of one dimensional, hypersmgular, mtegro-di referential 
equations of the second kind. As indicated previously, 
(47) must be augmented by a set of edge conditions if 
applicable. The authors are not aware of any theory that 
describes the solvability of (47). This subject is a mat- 
ter of ongoing research. However, for certain simple 
cases, one of which is described in the next section, (47) 
is greatly simplified and solvability theorems do exist. 

The characterization of the integral equation ker- 
nels by (48) greatly simplifies the numerical solution of 
(47). Analytical results for the logarithmic, Cauchy, 
and Hadamard kernels are available in many cases 11 and 
the continuous portions of the kernel can be computed 
by straightforward numerical integration. 

Results 

In this section, we consider the scattering of inci- 
dent sound by an infinitesimally thin cylindrical duct ot 
unit radius. The duct can be stationary or in motion. 
The interior and extenor walls of the duct are assumed 
to be hard. For this case, the complexity of the bound- 
ary integral equations is reduced significantly. 

Boundary Integral Equation Formulation 

From (12) and (18) the boundary conditions are 


" e "“’ r (lN )" (Z ' )dZ ’"° 


Z e[a, b]. (49) 


Differentiating (49) with respect to Z and using the re- 
lationships between total, scattered, and incident pres- 
sures yields the boundary equations 


Equation (51) is used below to simplify the field equa- 
tion (24). 

Define the jump in scattered pressure across the 
duct wail by 


AQ,(z) = Q;(z)-Q;(z) Ze[a,b]. 

(52) 

Referring to (25-26) obser/e that 


d-[f](r,Z) = -cT[f](r,Z) 

(53) 

and 


s-[f](r,Z) = s*[f](r,Z). 

(54) 


Thus, applying (51-54) to (24) produces the field equa- 
tion 

Q,(r,Z) = -d[AQ,](r.Z). (55) 

Therefore, the scattered acoustic pressure in the sound 
field is written as a double layer with density given by 
the scattered pressure jump. Since the interior and ex- 
terior duct surfaces are the same, we have omitted the 
superscript on the double layer operator. 

A single integral equation for AQ, is obtained from 
(49) as follows: Use (42) to evaluate the normal de- 
rivative of (55) on the exterior wall, then combine this 
result with the extenor boundary condition in (49) to 
give 


r —r 

j e M ' D,[AQ,](Z')dZ' = 



(Z')dZ' 


Z e[a,b] 


(56) 


It is advantageous to rewrite (56) as the system of 
equations 

D N [AQ,](Z) = ^| (z) Ze[a,b] (57a) 

and 
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g[AQ,] = C t> , 

(57b) 

C.-g^fy] 

a= !;<••■■<) r,i * 

where the functional g is defined by 


s '■ [i] 

where 

g[f]= 1 e M ’ D„[f](Z')dZ' 

-go 

(58) 

g ( ”- 0 [f] = g[(Z-a)' , (b-Z) i f] 


( 63 ) 


(64) 


and 


C 0 



- JQ 


(Z')dZ'. 


(59) 


Equation (57a) is obtained by differentiating (56), and 
(57b) by evaluating (56) at the trailing edge. No infor- 
mation is lost by this reformulation. For the stationary 
case, (57b) is satisfied trivially. By performing a local 
analysis on the kernel in (57a) it can be shown to have 
the form 


D n (Z-Z') = 

Bln|Z-Z'| + 


A 

(Z-Z') 2 + 

K.(Z-Z') 


(60) 


where A and B are known constants and the kernel K B 
is continuous and simple to evaluate numerically. 

To obtain a unique solution to (57a,b), the behavior 
of the pressure jump at the duct leading and trailing 
edges is required. It is known that the jump in pressure 
has the following asymptotic behavior: 


The notation in (64) will be used with other integral 
operators in the remainder of this analysis. Substituting 
(62-62) in (57a) gives the first kind integral equation 

K (/l - /l) [ Y ](z) = q(z) Z e[a, b], (65) 


where 


and 


K l/r/ - [y](z) = D, /r/ 2 '[y](z)- 


( 66 ) 


T 1 

q(Z) = 


fQi 

dr 


(Z)- 


L ' D ( f 2 , ' K) [l](z) Zs[a,b] 


(67) 


«^[1] 


The kernel for the operator in (66) has the same form as 
(60). After solving (65) for y , (63) is used to calculate 
a. 


AQ S = 


| o( VZ - a ) Z — > a 

i| 

V v b - 


+ Vb-Z Z -> b" 


■ (61) 


Based on (61) we assume a solution of the form 


Numerical Solution 

Due to the edge behavior associated with y , it is 
natural to expand y in a senes of Chebyshev polyno- 
mials of the second kind. Thus, constants (y • f are 
sought such that 


AQ,(Z)-C^ z _ a + 

1 

v ; (b-Z)(Z-a)y(z) 


(62) 


where a is an unknown constant and y is an unknown 
continuous function. Note that if the duct is stationary, 
then a = 0 , 

A method is now developed in which the determi- 
nation of y is separated from the calculation of a. 
This yields an integral equation for y that is relatively 
simnte to solve, Use (6 2) in i57b) to obtain 


y(z) = Zy,u J 

J-0 


2Z - a - b 


where 


Uj(x) = 


b-a 

sin [( j 1) cos' 1 x] 
sin(cos'' x) 


(63) 


(69) 


Golberg 9 has shown that if K is not an eigenfrequency 
of K, then a unique convergent expansion such as (68) 
exists for integral equations with kernels of the type 
given in (60). 
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( 75 ) 


To solve (65) numerically, we truncate the expan- 
sion in (68) and apply the collocation method. Other 
popular projection techniques, such as Galerkin’s 
method, require an additional numerical integration 
relative to collocation. With proper choice of colloca- 
tion points and numerical quadrature scheme for the 
continuous portion of the kernel, Golberg has shown 
that the accuracy obtained by Galerkin’s method and 
collocation are equivalent for this problem. Thus, col- 
location yields the same accuracy as Galerkin’s method, 
but with substantially less computational work. 

The numerical solution begins by choosing the 
number of terms in the expansion (68), N 0 +1. This 
number is a function of the modified wave number K. 
If 7 denotes the approximate solution, then 

•7tZ) = Lv^ (70) 

The collocation points, { Z j } ' are given by the ze- 
roes of the N 0 + l-th Chebyshev polynomial of the first 
kind and the numerical integration scheme tor the con- 
tinuous kernel is chosen as Gauss quadrature with 
weights and nodes based on second kind Chebyshev 
polynomials. Evaluating (65) at the collocation points 
yields the linear system 



[U,](Z J ) = q(z j ) 


j= +1 


(71) 


for the unknown expansion coefficients. The invertibil- 
ity of the linear system has been established by Golberg. 
To compute (71), integrals of the type 


i 

— Inix-x'ldx' xe(-co,l] 

V 1 — x ' 

-l 

are encountered. Analytical results are obtained for 
(72-75) by applying the Plemelj-Sokhotski theorem' and 
its logarithmic analog to the complex function 

R J (w) = [w-(w 2 -1 Y 1 

with branch cut 


| Re( w)J < 1 lm(w) = 0. 

The branch is defined such that (w“ — l) /fc is real and 
positive if w is real and greater than one. This analysis 
yields the following results: 


f Vl-x'^UjUOdx' 


(x-x')' 

Uj(x) xe[-l,l] 
(x + Vx : - l ) 


= -rt(j+0 : 


w* 1 




x<-l 


(76) 


|Vl-x' 2 U 0 (x')ln!x-x'|dx' = 

lT,(x)-In2 x«[-l,l] 

Z ji^ + Vx^T) - 

In ( — x t V x ” -l) x < -1 


(77) 


Vl- x' 2 Uj (x’)dx' 

(x-x') 2 


X 6 (-00,1], (72) 


j Vl-x' 2 U^xOlnix-x'jdx' xe(-oo,l], (73) 


k r 

J 


ltx' dx' 

* “ (x-x') 2 


x g(-oo, 1 ] , (74) 


and 


jVl-x^U^xOlnjx-x'idx' = 


1 j + 2 \ f ‘jV'v 

xe[-U] 

J + 2 j 


(x + Vx 2 -l) 

j + 2 

x < -1 j > 0 

J 
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Il + X' dx' 

V 1 “ x ' (x-x') : 

[0 x s[-l,l] 

d 1 X < — t 


(l-x)Vx 2 -1 


(79) 


and 


i/ir^ ln i x ~ x '! dx ' 


x + In 2 x <=[-!, l] 


x 4- In 2 - Vx : - 1 + 

In ( — x 4- Vx : - l ) 


x < -I 


( 80 ) 


In the above, ^ is the j-th order Chebyshev polynomial 
of the first kind 


In the first study, N B =20 point sources with tip 
Mach number = 1.2, were used to simulate the 
fan noise. Four forward flight Mach numbers (M F = 
0.0, 0.2, 0.4, and 0.6) were considered. Field calcula- 
tions in a plane perpendicular to the fan plane and paral- 
lel to the duct axis are presentee in figure 3a-d. 

An examination of the pressure fields inside the 
duct reveals that the wavelengths of axial modes propa- 
gating in the direction of motion decreases with increas- 
ing M F . A spectral analysis of the axial wave structure 
is beyond the scope of this work. However, the number 
of waves per unit length of the dominant axial mode 
present can be approximated by visual inspection. In 
table i, these observations are compared to the axial 
wave numbers for the first radial mode from classical 
spinning mode theory for ducts of infinite length. If 
v represents the smallest zero of the function J^, then 
the theoretical axial wave numbers are given by the 
formula 12 



Tj (x) = cos( jcos~‘ x). (31) 

Ex am ples 

The analysis leading to the BIEM is independent of 
the fan noise source description. For illustrative pur- 
poses, it is expedient to assume simplified source 
mechanisms with analytical expressions for the incident 
field. In the results presented here, a collection ofN B 
equally spaced point axial dipoles of unit strength lo- 
cated inside the duct at a radial distance 0.9r D and 
spinning with angular speed Q are used to simulate the 
loading noise produced by the fan (see figure 2). An 
analytical description exists for the incident field in- 
duced by this configuration 1 . The use of more sophisti- 
cated source processes is considered in the conclusion 
section. 

To demonstrate the versatility of the BIEM, several 
studies were conducted for the above problem. In each 
study, continuous portions of the total acoustic pressure 
field in the unstretched, moving frame were calculated. 
The numerical methods described in (66-81) were im- 
plemented on a Cray YMP computer at NASA Langley 
Research Center. For each set of parameters considered 
in the studies, the computational time for both field and 
integral equation calculations was 2-8 minutes. The 
acoustic fields displayed are composed of 20,000- 
20.000 ooserver point calculations. 


where is the m-th order Bessel function of the first 
kind. Only propagating modes are considered. Note 
that k L <0 corresponds to axial waves traveling from 
the fan face to the inlet. The computed resuits appear to 
be in agreement with theory. 

It is also noted that the angle between the line of 
peak noise and the duct axis decreases with increasing 
Mach number. This agrees qualitatively with the results 
of Rice, etal. 13 . 


m f 

K 

k,/2* 

(Theory) 

k*/2* 

(Observed) 

0.0 

24.0 

-1.44, + 1 .44 

-1.5, +1.5 

0.2 

24.5 

-2.47, +0.88 

-2.5 

0.4 

26.2 

-4.23, -0.59 

-4.0 

0.6 

50.0 

-7.59, -HL43 

-7.5, -0.5 


Table 1: Propagation Properties for Figure 3 Results 
(M-np = 1.2 L 0 =1.0 m = N 3 ) 

Using the same source configuration as above, the 
effects of increasing tip Mach number for fixed fright 
Mach number M F =0.8 were examined (see figure 4a- 
d and table 2). With regard to propagation and radiation 
characteristics, similar comments as in the previous 
study apply for cases 4a and 4b. 

In cases 4c-d, the waves moving forward in the 
duct are in agreement with the theoretical results. 
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However, waves traveling toward the exit are present 
that are not accounted for in the theory. The modified 
frequencies for these two cases are relatively close to 
some eigenfrequencies for the interior Dirichlet prob- 
lem. That is, the eigenifequencies are occurring at the 
zeroes of J^. There appear to be resonant radial modes 
present. Condition numbers for the linear system (71) 
at these frequencies increase significantly. The numeri- 
cal results are therefore questionable. This phenomenon 
has been well studied in the literature and is usually 
associated with ficticious interior eigenfrequencies 
while solving exterior problems. The interested reader 
is referred to the work of Kleinman and Roach 14 for a 
comprehensive theoretical discussion on the removal of 
the ficticious eigenfrequencies. In the present work, the 
intenor is real and the eigenfrequencies are not neces- 
sarily ficticious. Research into the subject is ongoing. 
The authors believe that the ill conditioning can be 
mitigated by the use of singular value decomposition 
methods. We further contend that acoustic treatment of 
the duct interior will eliminate the problem entirely. 


M-np 

K 

kj 2* 

(Theory) 

kjl-x 

(Observed) 

0.5 

16.7 

none 

none 

0.7 

23.3 

-3.06. -6,84 

-6.5 

0.9 

30.0 

-1.02,-11.71 ; 

- 1 1.5, *3.0 

1.1 

36.7 

-0.04, -15.52 

-15.0. +6.0 


Table 2: Propagation Properties for Figure 4 Results 
(M F = 0.8 L 0 = 0.5 m=N B ) 

The third example was chosen to demonstrate the 
capability of the BIEM to computer higher harmonics 
and different portions of the acoustic field (figure 5a-c). 
in this case, N a =lo, M F =0.2, and M tlp =L7. 
Sound pressure levels are plotted in an observer plane 
two duct diameters beneath the duct. The kinematic 
properties and the observer locations were selected to 
correspond to those used for tests conducted with the 
Langiey ducted propeller simulator 13 . Direct compan- 
ions with the results in reference 15 are not possible 
because of the simplified source model used for the 
BIEM. The results do show, however, the ability of the 
BIEM to compute the sound field in regions of interest. 


The results presented here demonstrate that the 
BIEM is a versatile and computationally efficient tool 
for predicting ducted fan engine noise. Qualitative ra- 
diation and duct propagation results can be obtained by 
using simplified source models such as spinning point 
or line sources. By tuning the strengths of the resulting 


monopoles and dipoles to account for fan loading and 
thickness effects, it is believed that the BIEM can pro- 
duce results that are useful for quantitative studies. 

Other boundary integral techniques have been de- 
veloped for the problem of scattering of incident sound 
by a thin duct with hard walls 10 . The BEEM developed 
here is valid for many situations of interest and features 
extensive mathematical analyses on the integral equa- 
tion kernels. The analyses yield expressions for singu- 
lar and logarithmic integral operators that can be evalu- 
ated in tenns of known functions and continuous por- 
tions that can be evaluated by simple numerical quadra- 
ture schemes. This versatility and depth ot analysis, 
absent in the referenced works, simplify the calculations 
considerably. 

Realistic duct geometry is included in the boundary 
integral equation formulation. Implementation of an 
arbitrary duct profile requires the solution of a system 
of two hypersinguiar integral equations. The inclusion 
of a duct centerbody produces another integral equation 
with the same properties. Both the duct profile and 
centerbody have interior regions that produce ficticious 
eigenfrequencies for the Neumann or Dirichlet bound- 
ary vaiue problems. This difficulty can be alleviated in 
several ways. The method of Burton and Miller 1 ap- 
pears to be adaptable to the BIEM presented here. Fic- 
ticious eigenfrequencies are not present if the engine 
components are acoustically lined. 

The use of passive noise control techniques are 
included in the BIEM. A spatially varying, locally re- 
active liner model appears in the boundary conditions. 
This property makes the BIEM attractive for active and 
passive noise control design studies. 


1. J. Lan: Acoustic Shielding from a Short Rigid Due: 
Around a Propeller, Master’s Thesis, The George 
Washington University, August 1993. 

2. J. Buhier: A New Boundary Integral Method tor 
Predicting the Acoustic Field Generated by Ducted Ro- 
tating Sources, Master’s Thesis, The George Washing- 
ton University, August 1994. 

3. R. Martinez: Aeroacoustic Diffraction and Dissipa- 
tion by a Short Propeller Cowl in Subsonic Flight. 
NASA CR 190801, April 1993. 

4. R. Martinez: A Boundary Integral Formulation tor 
Thm-walled Shapes of Revolution, J. Acous. Soc. Am., 
87(2), February 1990. 


I l 


American Institute of Aeronautics and Astronautics 




5. M.A. Hamdi and J.M. Ville: Development of a Sound Some Exterior Boundary Value Problems, Proc. Roy. 
? adizticn Model for a Finite Length Duct of Arbitrary See. Ser. A. 323, !9' ? L pgs. 201-210 

Shape, AIAA Journal, Vol. 20, No. 12, December 1982, 
pgs. 1687-1692. 

6. M.K. Myers: On the Acoustic Boundary Condition in 
the Presence of Flow, Journal of Sound and Vibration, 

71(3), 1980, pgs. 429-434. 

7. F.D. Gakhov: Boundary Value Problems, Dover 
Publications, Inc., New York, 1966. 

8. V.L Smirnov: A Course of Higher Mathematics, Vol. 

IV, Pergamon Press, Oxford, 1964. 

9. M. Golberg: The Convergence of Several Algorithms 
for Solving Integral Equations with Finite-Part Inte- 
grals, Journal of Integral Equations 5, 1983, pgs. 329- 
340. 

10. F. Farassat: Introduction to Generalized Functions 
with Applications in Aerodynamics and Aeroacoustics, 

NASA TP 3428, May 1994, (Corrected Copy - February 
1996) 

1 1. A.C. Kava and F. Erdogan: On the Solution of Inte- 
gral Equations with Strongly Singular Kernels, Quar- 
terly of Applied Mathematics, Volume XLV, Number 
1, April 1987, pgs. 105-122. 

12. Walter Eversman: Theoretical Models for Duct 
Acoustic Propagation and Radiation, Chapter 13, 

Aeroacoustics of Flight Vehicles: Theory and Practice, 

Vol. 2: Noise Control, NASA RP 1258, August 1991. 

13. E.J. Rice, M.F. Heidmann, and T.G. Sofnn: Modal 
Propagation Angles in a Circular Duct with Flow and 
their relation to Sound Radiation, 17th Aerospace Sci- 
ences Meeting Paper AIAA-79-0183, New Orleans, 

1979. 

14. R.E. Kleinman and G.E. Roach: Boundary Integral 
Equations for the Three-Dimensional Helmholtz Equa- 
tion, SIAM Rev., Vol. 16, No. 2, 1974. 

15. R.H. Thomas, C.H. Gerhoid, F. Farassat, O.L. Santa 
Maria, W.E. Nuckolls, and D.W. De Vilbiss : Far Field 
Noise of the 12 Inch Advanced Ducted Propeller Simu- 
lator, 33rd Aerospace Sciences Meeting and Exhibit 
Paper AIAA-95-0722, Reno, 1995. 

16. A.J. Burton and G.F. Miller: The Application of 
Integral Equation Methods to the Numerical Solution of 


American institute of Aeronautics and Astronautics 


12 



Figure 2: Geometry and Coordinate Definitions for 
BIEM Noise Prediction 
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Figure 3: Total Pressure Field 

m = 20 M ,, P = 1.2 L n = 1.0 
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Figure 4: Total Pressure Field 

m = 20 Mr = 0.8 L n =0.5 
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Acoustic Radiation and Propagation for an Acoustically Treated Engine Duct 
Via a Boundary Integral Equation Method 

M. H. Dunn and J. Tweed 
Old Dominion University 
Norfolk, VA 

F. Farassat 

NASA Langley Research Center 
Hampton, VA 

Abstract 

The prediction of ducted fan engine noise using a boundary integral equation method (BIEM) is 
presented. Governing equations for the BIEM are based on the assumption that the duct is approximated 
by a thin, finite length circular cylinder. Acoustic processes are assumed to be linear, generate spinning 
modes, and occur in a uniform flow field. A mixed boundary value problem (BVP) is defined that 
describes the scattering of incident sound by the cylinder. The boundary condition on the duct interior 
wall allows for an axially segmented locally reacting liner. Using potential theory, the BVP is recast as a 
system of hypersingular boundary integral equations for the unknown single and double layer potential 
densities. BIEM derivations and solution methodology are demonstrated for the scattering of incident 
sound generated by simple sources in a low speed uniform flow field. The key features of the BIEM are 
its computational speed and efficiency, versatility, validity over a wide range of frequencies, and the 
ability to compute rapidly any portion of the sound field without the need to compute the entire field 
Propagation results obtained are in agreement with the classical results of spinning mode theory for 
infinite ducts. Various calculations are presented to illustrate the utility of the BIEM as a tool for 
conducting active and passive noise control studies. 




Notation 


denotes that a quantity is dimensional when appearing over a variable 
a axial coordinate of duct trailing edge in stretched, moving frame 

b axial coordinate of duct leading edge in stretched, moving frame 

c,c ambient sound speed 

k — mM nr characteristic wave number of m-th circumferential mode 

L d ratio of duct length to duct diameter 

m circumferential mode number 

M F - V F /c forward flight Mach number 

M tip = r D Ojc tip Mach number 

N „ number of fan blades 

D 

p’ Eulerian description of total acoustic pressure field 

p’ Eulerian description of incident acoustic pressure field 

p' s Eulerian description of scattered acoustic pressure field 

(r, tp,z) cylindrical coordinates in stationary frame 
r D duct radius 

7 ,t time 

u' Eulerian description of acoustic velocity field 

V F , V F forward flight speed 

Z axial coordinate in stretched, moving frame 

o’ - specific surface acoustic admittance on duct exterior (interior) 

ft = yjl- M 2 f stretching parameter 




3 

k - kjP stretched characteristic wave number 

p 0 ambient density 

f2 shaft speed (radians/second) 

Introduction 

The reduction of tonal noise radiated by turbofans is a subject of ongoing aeroacoustics research. 
Advanced analytical tools for predicting the sound radiated from engine ducts can facilitate the design of 
active and passive noise abatement technology. To be useful in design studies, prediction tools should be 
fast, versatile, accurate, valid for a wide range of frequencies and engineering situations, and 
implementable on mainstream computer systems. The ability to compute any portion of the sound field 
without the need to calculate the entire field is an important attribute in this regard. Computational 
approaches such as Finite Element Methods, CFD, and Computational Aeroacoustics lack this property 
Therefore, their usefulness for parametric noise reduction studies is limited. 

In this paper, a Boundary Integral Equation Method (BIEM) for the prediction of ducted fan engine 
noise is presented. The work is an extension of the BIEM discussed in reference l***. The method is 
based on the equations of linearized acoustics with uniform inflow, and features analytical and 
computational techniques that minimize the consumption of computer resources. We validate the BIEM 
by reproducing qualitative propagation results from infinite duct theory. Various computational studies 
are presented to demonstrate the effectiveness of the BEEM as a tool in the study of active and passive 
noise control. 

It is assumed that the engine duct is approximated by an infinitesimally thin, circular cylinder of finite 
length. We adopt a scattering approach in which the acoustic pressure field is split into known incident 
and unknown scattered components. The acoustic source processes are assumed to generate an incident 
pressure field that can be represented by a superposition of spinning modes. 
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In a frame of reference moving with the engine duct and in regions of space not occupied by 
scattering surfaces, the modal components of scattered acoustic pressure satisfy Helmholtz’ equation. An 
exterior boundary value problem for the scattered pressure is obtained by including a far field radiation 
condition and boundary conditions on the interior and exterior cylinder surfaces. The boundary 
conditions allow for hard walls or a locally reacting liner that can be uniform or segmented 

We solve the boundary value problem by expressing the scattered pressure field as a sum of double 
and single layer Helmholtz potentials with unknown densities. Application of the boundary conditions to 
the layer representation yields a system of one-dimensional, singular boundary integral equations for the 
layer densities. For uniqueness, conditions on the densities at the duct leading and trailing edges must be 
imposed. The source terms for the system are related to the known incident pressure and its normal and 
tangential derivatives. The system of boundary integral equations, edge conditions, and method of 
solution comprise the BIEM. 

***CITE DIFFERENCES WITH PREVIOUS WORK ON INTEGRAL EQUATION METHODS 
FOR DUCTED FAN NOISE PREDICTION - MARTINEZ & THE FRENCH GUYS - NAMELY, 
DEPTH OF KERNEL ANALYSIS, SOLVABILITY RESULTS, NOTATION, DESIGN TOOL 

Boundary Value Problem Derivation 

We consider an engine fan surrounded by an infinitesimally thin, finite length circular cylinder 
translating in the +z (axial) direction with uniform speed V F . The fan is composed of N B equally 

spaced blades and rotates with shaft speed fl . The incident acoustic pressure field generated by the fan 
is given. Linear conditions are assumed to apply and the inflow is uniform. 

In the analysis that follows, all quantities are nondimensional: length by 7 D , mass by p,/ F , and time 

by Or' 





Governing Differential Equations 


5 


The total acoustic pressure in the sound field is split into known incident and unknown scattered parts 

p'(r,y/,z,t) = p;(r,(//,z,t) + p' 5 (r,y/,z,t). (1) 

In regions of space that contain no scattering surfaces, p\ is governed by the homogeneous wave 
equation 


" 1 & 1 d 

( d \ 

\ a 2 a 2 ' 

c 2 dt 2 r dr 

l (Jr ) 

1 r 2 cty 2 dz 2 


\p> o. 

Acoustic pressure and velocity are related through the acoustic momentum equation 

du' 


( 2 ) 


at 


+ Vp’ = 0. 


( 3 ) 


In a frame of reference moving with the duct, all dependent acoustic variables can be expressed as 
linear superpositions of spinning modes. For example, the scattered pressure has the form 

p'Xr, ¥ ,Z,t)= | (4a) 


and the total acoustic velocity is written 


u'{r,[f/,Z,t) = X(7 m (r.Zy 

m(N B =l 

where the stretched, moving axial coordinate Z is given by 

it. 




(4b) 


( 5 ) 


Incident and total acoustic pressures are written similarly. 

Modal amplitudes in the BIEM are calculated term by term. For notational convenience, the 
superscript m on the modal coefficients is dropped hereafter. 

Define the dependent variables Q, Q s , Q , and V r by 

Q(r,Z)= P(r,Z)e lM ' z (6a) 




QXr,Z)= P s {r,Z)e' M ' z 
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(6b) 


Q(r,Z) = ^(r,Z)e w 
V r (r,Z) = U r {r,Z)e‘ MFZ 

Combining (4-6) with (2) yields the two dimensional Helmholtz equation 


(6c) 

(6d) 


1 d 

( 

r dr 

J* 

dr) 


d 1 m 2 2 

H — + K~ 

dZ 2 r 


Qs= o 


(7) 


for the m-th scattered coefficient. Using (3-6a), the m-th radial component of the momentum equation 
(3) can be written as 


e' M ' Z V r (r,Z)-P M 


TIP 


M, 


F J 

— OO 


e‘ K,Z ^-(r,Z')dZ’. 


If the duct is stationary, then (7) and (8) reduce to 


\_d_ 

r dr 


< d} 

\ 


+ 


3 2 2 +K 

dz r 


P= 0 


and 




( 8 ) 


(9) 


( 10 ) 


Equations (7-10) are valid for points not lying on the stretched duct. 


Boundary Conditions 

To meet FAA noise certification levels, it is necessary to treat the engine duct with noise suppression 
devices. In this work, the duct treatment is modeled by a uniform or axially segmented locally reacting 
liner. We consider the case of a hard exterior surface and lined (or hard) interior. To simulate actual 
duct treatment, we impose hard wall conditions at and near the leading and trailing edges on the interior 


wall. 




Let f{r,Z) be an arbitrary field function. Define the surface function f'{Z) by 
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f ± (Z)= limf(r,Z) Ze{a,b). 

r-*l 


(ID 


In reference 2***Myers BC*** it is shown that if and a represents the piecewise constant specific 
acoustic surface admittance on the interior, then in the stretched, moving frame of reference, the 
boundary modal coefficients of velocity and pressure satisfy the boundary equations 

V;(Z) = 0 Zs(a,b) (12a) 

and 

iM F a{Z) d 


-K'(z) 


I -,- K -z 


+ - 


\t F 


Q 


fi 2 K dZ 

In the absence of flow, the boundary conditions reduce to 


(Z) = 0 Z &(a,b) . 


(12b) 



U;{z) = 0 z e(a,b). 

(13a) 

and 




-U;(z) + a{z)P {z) = 0 ze(a,b). 

(13b) 


We use the acoustic momentum equation (8) to eliminate the normal velocity from the (12): 


f 

\ cJr ) 


-e 


-wM SQ 


iM 2 F a(Z) d 2 


(Z) + - 

dr J V / P 3 M tip k dZ 2 


Z e(a,b). 

(14a) 

'"^ Z Q-\z) = 0 Ze{a,b ), 

(14b) 


and 


r < 
e 


' M ' Z ^-(/,Z')dZ' = 0 
dr 


-oo 


Note that (15) is satisfied trivially for M F -+0 and provides no information. 


(15) 


The thin duct approximation produces a surface with discontinuous tangents at the duct leading and 
trailing edges. Consequently, there are infinitely many solutions to (14-15). For uniqueness, we must 
constrain the acoustic pressure at the duct edges. At the trailing edge we impose the Kutta condition 
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lm[Q(Z)-Q(Z)] = 0 


(16) 


For M F >0 , the Kutta condition ensures infinite pressure at the leading edge. For physically reasonable 
solutions to exist we require the acoustic pressure to be integrable in any region of space about the 
leading edge. 

The boundary value problem is finalized by specifying the behavior of the acoustic pressure in the 
farfield. To ensure that the acoustic field consists only of outgoing waves at infinity we impose the 
Sommerfeld radiation condition 


Jim y[p\ 


dQ 

dp 


+ i«Q\ = 0. 


(17) 


Using (1) and assuming that the functions Q t and its derivatives are known, satisfy (17), and are 
continuous on the duct surface, we rewrite the boundary value problem in terms of the scattered pressure: 


\_d_ 

r dr 


d \ d 1 m 2 2 

r — + — 7 r + K 

dr ) dZ 2 r 2 




(18a) 


f dQ^ 

V dr ) 


(Z) = -^(/,Z) Ze(a,b) 


(18b) 


(dQ s \ {z) JdQ^ 


{ dr) 


v dr j 


(zb 


iM 2 F a(Z) 


P'kMtip 


■— 2 d 2 ( 

e M ' a 


dZ 2 


-i — Z 

e O. 


(Z) = 


iM^a{Z) j- s K Tf z &_ 
/PkM 


TIP 


dZ 2 


-X——Z 


e Q \l,Z) Ze(a,b) 


(18c) 


r k 


l,Z')dZ ' 


r lJL z ‘dQ t 


e Mf ^L(l,Z')dZ' 
dr 


(18d) 


lm[Q;(Z)-Q;(Z)] = 0 


Z-+a 


( 1 8e) 




J|0 s |<W<x 


9? any region in r - Z plane 


(180 
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Boundary Integral Equation Formulation 

In this section, we convert the boundary value problem (18a-g) to a system of boundary integral 
equations. 

Helmholtz Potential Representation 

The Green’s function for the two dimensional Helmholtz operator that satisfies the radiation condition 
( 1 8g) can be written as 


1 C e 

G(r,r',Z- Z') = — cosmi// dy 

2k J R 


(19) 


where 


R = yjr 2 +r' 2 -2rr’cosy/+{Z-Z'). (20) 

For an arbitrary surface function / , we define the single and double layer operators, S and d , by the 


equations 


and 


iftr,Z) = \f(Z’)ir,Z-Z')dZ' 


d[/](r,Z) = l/(Z')J(r.Z-Z')JZ’, 

a 

where the kernels s and d are given by 

s{ r,Z-Z') = G(r,l,Z-Z ') 


and 


d{r,Z-Z') = -^{r,l,Z-Z'). 


( 21 ) 

( 22 ) 

(23) 

(24) 


The integrals in (21-22) are well defined for points ( r,Z ) not on the stretched duct surface. 

We define additional field operators by calculating radial and axial derivatives of S and d Denote 
by S r the operator 

S ,[/]M) = 4:ifl r - z ) = ]/(z')s,(r.z-r)dZ' . (25) 

CT a 


where 




(26) 


s r (r,Z-Z’) = J^(r,Z-Z'). 

Define the field operators S z , S zz , d r , d z , d zz and associated kernels in a similar fashion 

Using results from Helmholtz potential theory (reference ***Courant & Hilbert***), the scattered 
pressure can be written as a sum of single and double layer Helmholtz potentials. 

QXr. z) = S[a \r, Z) + d[& ](r. Z) (27) 

Equations (18a) and ( 1 8g) are satisfied by (27). The layer densities Q, and 0 2 are unknown surface 
functions. Once the densities are determined, (27) is used to obtain the scattered field at any desired 
location. We will show later that Q, is related to the jump in the normal derivative of scattered pressure 
across the duct surface and Q 2 is related to the jump in scattered pressure. 

Surface Operator Notation 

In order to apply the boundary conditions (18b-c) to (27), it is necessary to evaluate directly the 
single and double layer potentials and their derivatives on the stretched duct surface. The resulting one 
dimensional surface operators have both singular and nonsingular parts. 

For Z e [a, b\ and sufficiently smooth / , we define the surface operators S and D by 


b 


S[/](Z) = ] f(Z')S(Z - Z')dZ' 

a 

(28) 

and 


D[/K z ) = ]/{Z')D{Z-Z')dZ\ 

a 

(29) 

where the kernels S and D are given by 


S{Z-Z') = s(J,Z-Z') 

(30) 

and 

D(Z-Z') = d{l,Z-Z’). 

(31) 


The operators S r , S z , S zz , D r , D z , D zz and associated kernels are similarly defined. 




Singular Kernel Analysis 
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All of the above kernels are singular for Z -Z’ - 0 . We list without proof the asymptotic behavior 
of the kernels for \Z - Z’\« I . The singular portions of the kernels are obtained by local analyses of the 

Green’s function (19) and its derivatives. All kernels are written as sums of singular and bounded terms. 
Bounded kernels are denoted by the superscript B . 


S(Z-Z') = -~ln\Z-Z'\ + S B {Z-Z') 

(32) 

S r {Z - Z') = ~j~ln\Z - Z'| + S B {Z - Z') 

(33) 

S z (Z-Z')= . 1 ,+S*(Z Z') 

zV ' 2k{Z’-Z) zV ' 

(34) 


I 4{K 2 -m 2 ) + l 

S 77 (Z-Z') = - + — - ln\Z - Z ' + S*(Z - Z’) 

2tv{Z - Z') 16k 11 


D(Z - Z') = - — /*jZ - Z'| + D b (Z - Z') 


/ 4(k 2 -m 2 \ + 3 , 

D(Z - Z') 7 + HZ - Z'\ + Df(Z - Z ') 

K 2k(Z-Z ') 16k 1 1 


D z (Z-Z') = — ; - + D z B (Z-Z’) 

z( ’ 4 k{Z'-Z ) zV ’ 


(35) 

(36) 

(37) 

(38) 


D zz {Z-Z') = 


l 4{K 2 -3m 2 } + 3 


4k(Z - Z') 2 


■ + 


32k 


/n)Z - Z'| + D zz (Z - Z') 


(39) 


The leading behavior for the kernels (32), (33), and (36) is logarithmic. Therefore, the associated 
operators are weakly singular. The leading terms for the kernels (34) and (38) are of the Cauchy type. 
While the kernels (35), (37), and (39) are of the strongly singular Hadamard type. Integrals with Cauchy 
and Hadamard kernels are divergent and must be interpreted in the finite part sense. 




13 

All of the above kernels have the logarithmic portions extracted. Integrations involving these terms 
are well defined theoretically but difficult to achieve numerically. This problem is mitigated by the 
development of analytical results for the associated operators. Examples of this, as well as analytical 
results for the Cauchy and Hadamard terms, are presented in a later section. Calculations involving the 
continuous portions of the kernels are performed by straightforward numerical integration. 

Laver Continuity Properties 

Using the above operator notation, we state well-known continuity properties for the single and 
double layers and their derivatives as a field point approaches the stretched duct surface from the interior 
or exterior of the duct. For sufficiently smooth /(Z) with Z s(a,b) , we have the following results (see 

reference ***Courant & Hilbert***): 


/ /ms[/](r,Z) = S[/](Z) 

r-»/ 

(40) 

lim&{f\r,Z) = +^/(Z) + D[/XZ) 

r->/ r 2 

(41) 

iimS,{flr,Z) = *!~ f (z) + S'[flZ) 

r~*l 2 

(42) 

towd r [/X r > Z ) = D r[/X Z ) 

r->/ 

(43) 

HmS z {/lr. z )=S r [/lZ) 

r-*r 

(44) 

tty *zz[/]( r > Z) = ^Zz[/K^) 

(45) 

limd z {/lr,Z)^^(Z) + D z [fJZ) 

(46) 

/™ d a [/Kr.Z) = + D^[/p) 

(47) 


By applying (40-43) to (27) we establish the previously mentioned relationships between the layer 


densities and scattered pressure surface quantities. 




q;(z)-q;(z) = -o 2 (z) 


14 

(48) 


dr 


M- 




(z) = -a(z) 


(49) 


Boundary Integral Equations 

Equations ( 1 8a-c) and ( 1 8g) are replaced by an equivalent system of integral equations for O l and 
Q 2 . Using the operator notation (28-29, etc ) and the continuity results (40-47), the boundary conditions 
( 1 8b-c) can be written as 


\ i+s r Ve,Kz)+D t [aXz) = z 


(50a) 


and 


<-e 




(3 3 kM tip dZ" 


< -,JL Z M 

■ w ' s 

n 


lQilz) + 


iM 2 F a(Z) d 2 


p 3 KM np dZ 2 


-! Z 


I + D 


lQ,lz) 


iMlajZ) ( ~ iJ ^ z A 
j 3 3 kM tip dZ 


(50b) 


<? Mr Q, V.Z) Zz{a,b) 


where I is the identity operator. 

To obtain a unique solution of the system (50a-b), we must establish requirements on the layer 
densities at the duct leading and trailing edges. This information is provided by the as yet unsatisfied 
BVP equations (18d-f). Applying the layer representation for the scattered pressure to (18d-f) yields 


a 



QM) = 0 


(50d) 


and 


J|s[GXr.Z) + d[&](r.Z)|<»! < 00 


9? any region in r - Z plane . 


(50e) 
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The boundary value problem (18a-g) is thus replaced by the equivalent boundary integral equation 
formulation (50a-e). Equations (50a-b) comprise a system of strongly singular integro-differential 
equations 

Elaborate techniques are required to solve (50). For clarity, we consider the solution of (50) under 
simplified circumstances of practical interest. For a = 0 (hard interior wall) in (50b), we obtain the 
result Qi(Z) = 0 . This yields the classical Neumann boundary integral equation 


D,[ Q,\Z) = -^ r (l,Z) 2 € (a. A) 


(51a) 


and auxiliary conditions 


<3 



(51b) 


QA a ) = 0 . 


J|d[&X'. Z)|t/9?<ao 


9? any region in r — Z plane . 


(51c) 
(5 Id) 


Note that (51b) is satisfied trivially for M F = 0 . 

The solution of (51) for cylindrical ducts has been studied in the literature recently (references ***- 
***). The referenced works differ in their derivations, depth of kernel analyses, and solution techniques. 
In reference ***AIAA paper***, the analyses and techniques similar to those presented here were 
applied to (51). For the first time using boundary integral techniques, continuous portions of the acoustic 
field were calculated quickly (several minutes on a PC) and accurately. This was made possible by 
applying advanced analytical and numerical methods to the singular integrals in (5 1). 

Also of interest are “small” inflow Mach number situations. Such conditions occur at take-off and 
landing for example. Ignoring terms in (50) smaller than 0{l) for M F —*0, yields the approximate 


system 
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-j I+s r )[a](z)+D,[ap)=-^(/.z) z €(<.,*) 

| + ^s([ 0 iz)-^(|, + D) fe Kz)^ a(/ . z) 

***CHECK THIS ONE OUT FOR THE HARD WALL CASE*** 


Z &{a,b) 


(52a) 

(52b) 


S r[aX'.a-) + d r [a](/,a ) = -^(/.a) . (52c) 

&(«) = 0 • (52d) 

JW&X r,Z) + d[Qj](r,Z)|</SR <<» 9 ? any region in r- Z plane . (52e) 

9 ? 


In (52c) the notation a implies the limit as Z approaches a from the left. 

In the next section, the solution of (52) will be examined in detail. The validity of the approximate 
equations for 0 < M F < 0.4 *** THIS RANGE MAY CHANGE - NEED SOME HARD WALL 

CALCULATIONS TO VERIFY*** is demonstrated in the results section. The solution of the full 
equations (50) will be considered in future research by the authors. 
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List of Symbols 


a axial coordinate of duct trailing edge in moving frame 

b axial coordinate of duct leading edge in moving frame 

c ambient sound speed 

k = mN a M T]P nondimensional characteristic wave number of m-th circumferential mode 

L d ratio of duct length to diameter 

m circumferential mode number 

M = Vjc flight Mach number 

M np - r D Q/c tip Mach number (based on duct radius) 

N B number of fan blades 

N L number of liner segments 

N obs number of observers for TBIEM3 D output 

P total acoustic pressure 

P m m-th circumferential coefficient of total pressure 

P t incident acoustic pressure 

P m m-th circumferential coefficient of incident pressure 

P scattered acoustic pressure 

P m m -th circumferential coefficient of scattered pressure 

(r, y/, Z) cylindrical coordinates in frame of reference attached to duct 
r n duct radius 

r„ radial coordinate of spinning point dipoles 

i 


time 
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r 



a(Z) 



P 

K 


Po 

ct(Z) 

4(Z) 

Q 


t hoist from fan 
duct speed 

axial locations of liner segments 

= p i) c\^{Z)-i<y{Z)\ segmented, specific acoustic admittance on interior duct wall 
piecewise specific acoustic admittances 

= •>//- M : compressibility (stretching) parameter 
= kj p nondimensional stretched characteristic wave number 
ambient density 

segmented, acoustic susceptance on interior duct wall 
segmented, acoustic conductance on interior duct wall 
shaft speed (radians/second) 



Introduction 


This document describes the ducted fan noise prediction computer program TB1EM3D (Thin duct. 
Boundary Integral Equation Method, 3 Dimensional) The scattering of fan generated noise by a finite 
length, infinitesimally thin circular cylinder in a uniform flow field is considered The program, based on a 
boundary integral equation method (BIEM), calculates circumferential modal coefficients of the acoustic 
pressure at user specified field locations TBIEM3D features include versatility, rapid calculations, and 
ease of use. Theoretical and computational details can be found in references 1-4 

In a frame of reference attached to the duct, the fan generates spinning acoustic modes The thrust 
component of fan loading noise is approximated by a collection of spinning point thrust dipoles. A 
precise mathematical representation for the acoustic field due to this configuration has been implemented 
In many cases, TBIEM3D can be easily adapted to accommodate other source fields 

TBIEM3D employs cylindrical coordinates in a frame of reference attached to the engine (figures 1- 
3). The coordinate origin is at the center of the fan disc. The fan and duct translate in the +Z (axial) 
direction with uniform speed V . N B equally spaced blades comprise the fan. The shaft rotates with 


speed n (figure 2). 

The total acoustic pressure in the sound field is split into known incident and unknown scattered 
parts: 

P’{r,y/,Z,t) = P’{r,y/,Z,t) + P/(r,y/,Z,t) (1) 

Assuming linear conditions, all dependent acoustic variables can be expressed as superpositions of 

spinning modes For example, the scattered pressure has the form 


P;{r.ys,Z,t)= X/r(r,Z)e"" v ' (Q '^ . (2) 

m-- -a- 

Incident and total acoustic pressures are written similarly Modal amplitudes are calculated term by term 
The TBIEM3D code must be run separately for each desired mode 



The duct exterior is hard and the interior may be hard or lined. Passive noise treatment is modeled by 
an axisymmetric, locally reactive, segmented liner with user specified admittances The definition of 
specific acoustic admittance used by TBIEM3D 

a(Z) = p„c{4(Z)-io(Z)} (3) 

is consistent with the time factor in (2). Regions of the duct interior near the leading and trailing 

edges are assumed hard (figure 3) Any interior wall segment may also be rigid. These comments are 
summarized by the equation 

0 Z eja,Z,]^[Z v ,,b\ 

a(Z) = \ , \ < 4 ) 

[a, Zg(Z ] ,Z^,) j = l,...N L 

BIEM methodology is a three step process: Step 1) A Helmholtz boundary value problem (BVP) for 
the modal coefficients in (2) is derived. Step 2) Using layered Helmholtz potentials, the BVP is 
converted to a boundary integral equation formulation that features a set of hypersingular integral 
equations for the unknown Helmholtz layers. Step 3) The integral equations are solved and the acoustic 
field calculated from the Helmholtz potential representation. 

The TBIEM3D code is written in the FORTRAN programming language and employs IMSL 
mathematical library routines. TBIEM3D should be implementable on any computer that can 
accommodate FORTRAN and IMSL Some code modification may be required. For minimally adequate 
computational performance, a Pentium 133 processor (or equivalent) with 32 megabytes of RAM is 
recommended 

TBIEM3D input is relatively simple. Geometric, kinematic, and liner parameters are required. If a 
source description other than the one described above is desired, then the user must supply FORTRAN 
subroutines for the calculation of the incident field and its radial derivative Output from TBIEM3D 
consists of the modal coefficients of the complex pressure components [see equations (1-2)] at user 
specified field points Postprocessing of results is left to the user 
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The key feature of TBIEM3D is the ability to compute any portion of the sound field without the 
need to calculate the entire field. Competing methods such as finite differences and finite elements lack 
this property Other positive attributes include reduced consumption of computational resources, 
enhanced numerical accuracy, versatility, coupling of radiation and propagation both forward and aft, and 
validity over a wide range of frequencies. Consequently, the TBIEM3D code is well suited for 
parametric calculations. Many engineering studies of interest can be handled by TBIEM3D 

Questions, comments, and requests for discussions should be addressed to mhd3l4®aol.com 



Limitations and Comments 


1) At present, the TBrEM3D code can treat ‘'small” Mach number inflow Results obtained for 
A/ > 0.4 may be questionable TBIEM3D with no inflow restrictions will be made available when 
complete 

2) For large values of k , TBIEM3D computational time and storage requirements can increase 
considerably. Therefore, at typical fan operating conditions, it is recommended that the user calculate 
a maximum of three circumferential modes. Efforts are underway to improve TBIEM3D 
performance for high frequencies. 

3) It is well known from the theory of wave propagation in an infinite, hard walled duct that resonance 
occurs at certain discrete frequencies At these eigenfrequencies, the infinite duct problem is 
unsolvable. Theoretically, the finite, hard walled interior duct is solvable at all frequencies. Ill- 
conditioning in the TBIEM3D numerical system, however, is experienced at and near the infinite duct 
eigenfrequencies. TBIEM3D results at these eigenfrequencies show evidence of resonance but 
appear plausible. The numerical correctness of TBIEM3D at resonance has not been established. 
Therefore users should examine TBIEM3D results carefully when the hard wall interior option is 
activated 

4) For some applications, it may be convenient to place the sources outside the duct This is easily 
achieved with TBIEM3D The user must have either a > 0 , b < 0 , and/or r 0 > r D . 

5) Since the duct is approximated by an infinitely thin cylinder, the acoustic pressure is discontinuous 
across the duct surface Consequently, evaluation of the acoustic pressure on the duct wall is 
ambiguous It is recommended that if the pressure on the interior duct surface is required, then the 
user should place the observer a small distance off the duct toward the interior 
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Operating Instructions 


TBIEM3D operating parameters consist of a one line identifier, output file name and path, and 
physical parameters The code generates one output file containing values of program parameters and the 
complex modal coefficients of incident, scattered, and total pressure at user specified field points The 
output file is associated with logical unit 9 Access of unit 9 elsewhere in the calling program can lead to 
errors and should be avoided COMMON statements in TBIEM3D should be examined to avoid 
conflicts with the user program 

To activate TBIEM3D, the user’s calling program must have the FORTRAN statement 


CALL TBIEM3D( idem , outfile , m, N a , RPM , r D , a , b , r 0 , T, c, p 0 , V 


( 5 ) 


. rc • w:, ■ "o.s.vx,: .[z,],: > 


The notation |x, j' in (5) denotes a one dimensional array of length N Variables in the argument list 

are described in table 1 SI units are required for dimensional variables. 

The TBIEM3D output file contains the case identifier and program parameters followed by N OBS 
formatted lines containing the dimensional (pascals) complex modal coefficients of incident, scattered, 
and total acoustic pressure For each observer point, TBEEM3D writes the pressure components 
according to the following FORTRAN statements: 


WRITE(6,600) r , Z , Re[P , m ] . Im[P , m ] , Re[P s m ] , In{P ; m ] , Re[P m ] , Int[P m ] 

600 FORJVLVT(8EI 1.4) . 
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Examples 


Three examples are presented in this section to demonstrate TBIEM3D features and usage 
Kinematic parameters were chosen to simulate actual ducted fan engine operating conditions 
( \'i F = 0.40 , M T jp = 1.22, N 8 = 20, L D = 0.50 ) The examples differ in the acoustic treatment on 

the duct interior. Admittances for the three cases are given below. All calculations were performed on a 
Pentium 133 laptop computer with 32 megabytes of RAM. Graphical results displayed here are not part 


of TBIEM3D. 

For each of the three cases, two dimensional portions of the sound field are computed. Acoustic 
pressure and sound pressure level contours for the first modal coefficient are plotted in figures 5-7 The 
specific acoustic admittances used for the calculations do not necessarily correspond to actual conditions 
but were chosen for demonstrative purposes. Figure 4 contains the FORTRAN calling program that 
generated the results for figures 5-7 

Example 1 Hard inlet and hard exhaust. See figures 5a-b. Three minutes computational time required 
for 1 5000 field points 
Admittance: a(Z) = 0 Z e[a,Z>] 

Example 2 Lined inlet (one segment) and hard exhaust. See figures 6a-b. Five minutes computational 
time required for 15000 field points 


Admittance: a(Z) = 


[/-/ Z <=(0,0.475) 
1 0 elsewhere 


Example 3 Lined inlet (one segment) and lined exhaust (one segment). See figures 7a-b Five minutes 
computational time required for 1 5000 field points. 


Admittance a(Z) - 


1-i Z e(0, 0.475) 
0.5 Z (=(-0.475,0) 
0 elsewhere 
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In figures 8a-b, comparisons between the three cases in both the nearfield and farfield are displayed 
Sound pressure levels for the first modal coefficient on an arc of 200 field points about the duct center 
are calculated The radius of the arc is ten meters for the farfield example (figure 8a) and one meter for 
the nearfield example (figure 8b) and extends from the forward duct axis to the aft duct axis The results 
are plotted to show the effects of passive noise treatment. Calculations required approximately one 
minute. 
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Variable 

Description 

Comments 

ident 

Character* 80 variable, case identifier 

80 ASCII characters maximum 

outfile 

Character*80 variable, output file name 

80 ASCII characters maximum, file path 
may be included 

m 

Integer, Circumferential mode number 

See (2) and limitation 2) 

Nb 

Integer; Number of fan blades 

See figures 1-3; N B > 0 

RPM 

Real, shaft speed (revolutions per minute) 

Q = (nRPM)/30 

r D 

Real; duct radius (meters) 

See figures 1 -3 

a 

Real; axial coordinate of duct trailing 
edge (meters) 

See figures 1 -2 

b 

Real; axial coordinate of duct leading 
edge (meters) 

See figures 1 -2 

r o 

Real; radial location of spinning dipoles 
(meters) 

See figures 2-3 

T 

Real; thrust from fan (kilonewtons) 


c 

Real; Ambient sound speed (meters per 
second) 


Po 

Real; Ambient density (kilograms per 
cubic meter) 


V 

Real; Engine speed (meters per second) 

Vjc < 1 ; see limitation 1 ) 

V L 

Integer, number of liner segments 

For hard wall interior set N L =0 in 
which case limitation 3) may apply 

w 

Real; axial coordinates of liner segments 
(meters) 

See (3-4) and figure 3; if N. = 0 , then 
omit; a < Z < Z^, <b i = J....N L 

w;- 

Complex; acoustic admittances (non 
dimensional) for segmented liner 

See (3-4) and figure 3; if A r = 0 , then 
omit; some segments may be hard, i.e., 
<*, =° 

iV 

iy OBS 

Integer; number of observers for output 

0<N obs ^10 5 , Large values of N cm 
can lead to excessive computational time 


Real, radial coordinates of observer points 
(meters) 

See limitation 5) 

{z,Yr 

Real; axial coordinates of observer points 
(meters) 
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Figure 1 : Duct Geometry and Coordinate Definitions 
Cylindrical Frame Fixed to Duct 






Duct Exterior Wall 



Figure 3: Duct and Liner Geometry (Side View) 
Two Piece Segmented Liner 


Figure 4: Calling Program for TBIEM3D Examples 1-3 


c TBIEM3D examples for user document 


program myBIEM 

dimension 20BS (100000) , rOBS (100000) 

dimension Z liner (100) 

complex admit (100) 

character* 8 0 ident, out file 

ident = ' TBIEM3D Example #1* 

outfile = ' c:\myBIEM.txt' 

pref - 2. *10.** (-5. ) 

nharm = 1 

nblades = 20 

rpm = 3500. 

radiusd = 1.0 

zte = -0.5 

zle = 0.5 

radiusO = 0.95 

thrustc = 27.0 

sposnd = 300. 

density = 0.4 

V = 120. 


c 

c 

c Define observer points for r-Z field 
c calculations 

c 

c 

nz = 150 
nr = 100 
NOBS = nz*nr 
zinin = -2. 
zmax = 2 . 

dz = (zmax-zmin) /float (nz-1) 
rmin = 0 . 
rmax = 3 . 

dr = ( rmax-rmin) /float (nr-1) 
kount = 0 
do i = l,nz 

z = zmin+ (i-1 . ) *dz 
do j = l,nr 

kount = kount+1 
r = rmint ( j-1 . ) *dr 
ZOBS (kount) = z 
rOBS (kount) = r 
enddo 
enddo 


c 

c 

c 

c 

c 


Example 1: Hard Wall Interior 


NL = 0 

call TBIEM3D ( ident, outfile, nharm, nblades, rpm, 

1 radiusd, zte, zle, radiusO, thrustc, sposnd, 

2 density, V, NL, Zliner, admit ,N0BS, ZOBS, rOBS) 



Figure 4 (Continued): Calling Program for TBIEM3D Examples 1-3 


c Begin postprocessing for user document 
c graphics - NOT part of TBIEM3D 
£** + + + + + + + ** + + * + * + * + + •*■ + + + + + + + + + + + + + ★★ + ★ + + + ★★ + 
open (unit =10, file=' c: \ fort 10. txt 1 , status= 

1 * unknown 1 ) 

write (10, *) 'zone t = "1", i = f ,nr, , / j = ',nz 
rewind (9) 
do j = 1,15 
read ( 9, * ) 
enddo 

do iobs = 1 , NOBS 

read (9, 600 ) zz , rr , res , ais , rei, aii, ret, ait 

pmag = sqrt ( ret**2+ait**2+l . e-12 ) 

spl = 10 . *alogl0 (pmag/pref ) 

write { 10, * ) zz, rr, spl, ret 

enddo 


c * + **** + + + + + * + *** + + + + *** + ****>** + + + + + + * + + + * + + + + + + ** + 
c End postprocessing 

^'*- + + ****”*"*”* + -*- + * + *'* , *** + '*’'*' + *^‘*"*‘*‘* , + '*’' A ’ + '* , '*''*"*'****'*'*' A "*"*''*'** 
c 

c 

c Example 2: Lined inlet with one segment 

c 

c 

NL = 1 

Zliner ( 1 ) =0 
Zliner ( 2 ) = 0.475 
admit (1) = cmplx ( 1 . , -1 . ) 
rewind ( 9 ) 

call TBIEM3D (ident, out file, nharm, nblades, rpm, radiusd, 

1 zte, zle, radius 0, thrust c, sposnd, density, V, NL, 

2 Zliner , admit, NOBS, ZOBS, rOBS) 

c *** + ★ + + **★*★****•*• + + + ** *★★*★*** + *★★★★★*★★*★ + ***★★★* + 
c Begin postprocessing for user document 
c graphics - NOT part of TBIEM3D 
^ + + + + + + + + ■*' + + + + + + + + + + + + + + + + + + + + + + + + * + + + + + + + + +' + '*'***** * 
write <10, *) 'zone t = "1", i = '^r,', j = f ,nz 
rewind ( 9 ) 
do j = 1,15 
read ( 9 , * } 
enddo 

do iobs = l,NOBS 

read ( 9 , 600 }zz,rr,res,ais,rei,aii,ret,ait 

pmag = sqrt ( ret**2+ait**2+l . e-12 ) 

spl = 10 . +alogl0 (pmag/pref ) 

write ( 10, * ) zz, rr, spl, ret 

enddo 

^ + + + + + + + + + + + + + + + + + + + + + + + + + +■ + + + + + + + + + + + + + + + + + + * + + * + + + + + 
c End postprocessing 

c + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + * + + + + + + + + * + + + + + + + * 

c 

c 

c Example 3: Lined exhaust with one segment and 



Figure 4 (Continued): Calling Program for TBIEM3D Examples 1-3 


c lined inlet with one segment 

c 

c 

NL = 2 

Zliner(l) = -0.475 
Zliner{2) - 0. 

Zliner(3) = 0.475 
admit (1) = cmplx (.5,0.) 
admit(2) = cmplx ( 1 ., -1 . ) 
rewind { 9 ) 

call TBIEM3D { ident , out file, nharm, nblades, rpm, radiusd, 

1 zte, zle, radius 0 , thrustc, sposnd, density, V, NL, 

2 Zliner, admit , NOBS, ZOB5, rOBS) 

q+ + ** + + -A-******** + ** + '* _ + ***********'*'** > *‘* , *'* , ** , *'* , **'****'**'* , *‘* 

c Begin postprocessing for user document 
c graphics - NOT part of TBIEM3D 

+ + + + + + * + ★ + + * + * + + * + * + + 

write ( 10 ,*)* zone t = "l", i = 1 /nr, 1 , j = 1 , nz 
rewind ( 9 ) 
do j = 1,15 
read ( 9, * ) 
enddo 

do iobs - l,NOBS 

read {9,600}zz,rr,res,ais, rei , aii, ret, ait 

pmag — sqrt ( ret * *2+ait* *2+1 . e- 12 ) 

spl = 10 . *alogl0 (pmag/pref ) 

write ( 10, * ) zz, rr , spl, ret 

enddo 

£+ + ** + + + **** + * + * + + + + + + -*-**' + * + * + , *’*'*'* + + '*' 
c End postprocessing 

stop 

600 f ormat ( 8el 1 . 4 ) 

end 




Figure 5a: BIEM Example #1 - Hard Inlet and Hard Exhaust 

Acoustic Pressure 



Sound Pressure Level, dB (Re 20 |x Pa) 



Sound Pressure Levels 
0.40 M 11p =1.22 L d = 0 



Inlet and Hard Exhaust 



Sound Pressure Level, dB (Re 20 \ x Pa) 



Figure 6b: BIEM Example #2 - Lined Inlet and Hard Exhaust 
Sound Pressure Levels 

M c = 0.40 M-- = 1.22 L„ = 0.5 N = 20 



Re[p] (Pascals) 



Figure 7a: BIEM Example #3 - Lined Inlet and Lined Exhaust 

Acoustic Pressure 



Sound Pressure Level, dB (Re 20 p Pa) 



Figure 7b: BIEM Example #3 - Lined Inlet and Lined Exhaust 

Sound Pressure Levels 

M c = 0.40 M_ d = 1.22 U = 0.5 N D = 20 



Hard Inlet and Hard Exhaust 
Lined Inlet and Hard Exhaust 



Figure 8a: Effect of Passive Noise Treatment 
on Farfield Radiated Sound 
D.4 M T1D =1.22 N d = 20 L n = 0.5 r™ Q = 10m 



Hard Inlet and Hard Exhaust 
Lined Inlet and Hard Exhaust 
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Figure 8b: Effect of Passive Noise Treatment 
on Nearfield Radiated Sound 
0.4 M tid = 1 .22 N q = 20 L n = 0.5 r^ 




A Study of Noise Radiation from a Two 
Dimensional Scarf Inlet 

M.H. Dunn and Richard St. John 
Old Dominion University, Norfolk, Virginia 


1.0 Introduction 

The scarf inlet and exhaust have been recommended for the control of the 
noise radiation directivity from a ducted fan engine. There has been some 
experimental evidence that a scarf inlet does affect the directivity of the noise 
radiation pattern from the inlet and therefore, it is a promising passive noise 
control tool for aeroengines. Recently a computer code was developed at the 
Old Dominion University based on the boundary integral equation method 
(BIEM) for noise radiation from finite thin two dimensional ducts which 
allows the study scarf inlets and exhausts. The noise source inside the duct is 
taken as a point dipole source. We study the effect of the scarf inlet on noise 
radiation pattern from the inlet and the exhaust of a short duct with a conven- 
tional exhaust. The duct wall is assumed rigid and only one configuration of 
scarf inlet is studied. It is shown that the directivity of noise radiation from 
the inlet and exhaust is influenced by the presence of scarf inlet over a large 
range of frequency. The main part of the present study is for a stationary duct. 
However, we present some results at duct forward Mach number of M = 0.4 
which indicate that the scarf inlet also influences the noise directivity even 
when the duct is in motion. 

2.0 The Scarf Inlet 

Figure 1 shows the two dimensional scarf inlet used in this study. Note that 
the upper wall is longer than the lower wall. Only one configuration of scarf 
inlet with conventional exhaust was studied. The design is based on the long- 
est and the shortest axial nacelle length of the Langley 12 inch ducted fan 
model used in the scarf inlet test. This figure is to scale and the three circles 
of radii 0.2, 2.0 and 6.0 (nondimensionalized with respect to the duct width) 
indicate the arcs on which the acoustic pressure levels were calculated. The 
dipole point source is at the center of these circles and is located at the axial 


A Study of Noise Radiation from a l wo Dimensional Scart lnletluly 1, 1907 


1 




position of the rotor center in the model test. The dipole source models the 
thrust of the rotor with the dipole axis along the duct axis. The axial position 
of the dipole corresponds to the rotor position in the Langley 12 inch ducted 
fan model. 

3.0 The Acoustic Code (TBIEM2D) 

The acoustic code used in this study was developed at the Old Dominion Uni- 
versity. It is called TBIEM2D (Thin-duct BIEM 2 Dimensional). The duct is 
assumed finite in length, infinitely thin and can be in uniform forward 
motion. The interior of the duct can be lined and the two duct walls do not 
need to be of equal length. The liner can be segmented. In the examples given 
here we have used rigid duct wall option. 

The duct propagation and radiation problem is modeled as follows. The 
acoustic pressure is written as the sum of the incident and scattered pressures. 
The scattered acoustic pressure in the frame fixed to the moving duct satisfies 
the convected wave equation which is Fourier transformed in time into a 
Helmholtz partial differential equation. The boundary condition for this 
equation is obtained from the incident wave and the locally reacting liner 
property. The problem is then reformulated as a coupled system of one 
dimensional hypersingular linear integral equations with single and double 
layer sources of unknown strengths on the duct wall. The radiation condition 
is automatically satisfied in the derivation of these equations but additional 
conditions are imposed to insure uniqueness and physical validity of the solu- 
tion. 

The system of the integral equations are solved by the collocation method for 
the unknown source strengths of the single and the double layer sources on 
the duct wall. The scattered acoustic pressure is then calculated at any point 
in the radiation field from integrals on the duct surface whose integrands 
depend on the known strengths of the single and the double layer sources. 
This is essentially equivalent to using the Kirchhoff formula for moving sur- 
faces. The method of finding the scattered acoustic pressure is known as the 
boundary integral equation method (BIEM). The scattered pressure is added 
to the incident pressure to get the total acoustic pressure at the observer point. 

One of the most important advantages of our approach is that the acoustic 
pressure can be calculated at any point in the field without the need for calcu- 
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lating the acoustic pressure elsewhere. This is in contrast to the finite differ- 
ence and finite element methods which require the gridding of the space 
between the duct and the observer position and calculating the acoustic pres- 
sure at each grid point as well as at the observer position of interest. Thus the 
boundary integral equation method results in a very efficient computation 
tool for the acoustic pressure. The code TBIEM2D is run on a PC and has a 
simple input for point source excitation. The entire acoustic radiation field 
can be computed at a fixed observer time or at many observer times and then 
animated to see the propagation of the acoustic waves in the medium as well 
as the region of intense acoustic intensity and shadows. 

4.0 Scarf Inlet Acoustic Results 

We present our results below using nondimensional length and wave numbers 
as follows. All lengths are nondimensionalized by the duct width 

W = 0.30 m. The nondimensional wave number is defined as k — c oW/ c 
where co and c are the radian frequency and the speed of sound taken as 340 
m/s. In all the calculations presented, here, the source strength of the dipole 
inside the duct is kept fixed. Therefore, the relative levels of the acoustic pres- 
sure for the scarfed and conventional (unscarfed) inlet remain unchanged at 
the same nondimensional distance from the duct. 

Figure 2 presents the directivity of the acoustic pressure level at the distances 
0.2, 2.0 and 6.0 from the source and at the wave numbers k = 18 for both the 
scarf and conventional inlets at M p — 0 . Figures 3 and 4 show similar 

results at k = 27 and 36, respectively. We first consider the results of Fig. 2 at 
k = 18. The circle at the radius R = 0.2 is inside the duct and the directivity of 
the acoustic pressure level is clearly very similar to an axial dipole in the free 
field for both scarfed and unscarfed inlets. In the inlet radiation field, the scarf 
inlet radiation pattern becomes asymmetric with respect to the duct axis with 
a slight noise reduction in the upper and lower quadrants at angles of about 
60° at radii R = 2 and R = 6. At R = 6, there is a slight increase in the acous- 
tic pressure level of the noise radiated from the exhaust in the lower quadrant 
from the duct in scarf inlet. At k = 27 and k = 36, there are more lobes in the 
inlet and the exhaust radiation patterns of the two ducts. The asymmetry of 
the radiation pattern of the scarf inlet, as compared to the conventional inlet, 
is clearly visible in Figs. 3 and 4. In the inlet radiation field, there is a reduc- 
tion of the noise levels in the upper quadrant and an increase in the lower 


A Study of Noi.se Radiation from a Two Dimensional Sear! InletJuly 1. 19 ( )7 


3 




quadrant for the scarf inlet as compared to the unscarfed inlet. This effect is 
seen at large angles to the duct axis. There is some small changes in the 
exhaust radiation field by the scarf inlet at large angles to the duct axis. 

Figures 5 to 7 show the directivity of the acoustic pressure level at three 
radial distances from the source and at three wave numbers for the scarf and 
conventional inlets at M p = 0.4. The effect of the motion of the duct on the 

radiation pattern inside the duct is now much more pronounced than at 
Mp - 0. In the exterior of the duct and for the inlet radiation, there is a 

reduction on of the noise level in the upper quadrant and an increase in the 
noise level in the lower quadrant for the scarf inlet as compared to the con- 
ventional inlet at large angles from the duct axis. More surprising is the 
change in radiation pattern from the exhaust for the two inlet designs, partic- 
ularly at k = 18. This may have important implication for the a real engine 
noise radiation. 

As suggested by Dr. Robert Dougherty of the Boeing Commercial Airplane 
Group, a study of the relative acoustic power radiated in different quadrants 
can be useful in assessing the effect of scarf inlet on the control of noise radi- 
ation from ducted fans. We define an indicator of the acoustic power by the 

following relation for M p = 0 : 


P = 


fVo + rc/2 

Vo 


p ,2 (R, 


where \y () is 0, n/2 , n and 3tc/2 for the first to fourth quadrants, respec- 
tively. Note that the first and the fourth quadrants correspond to the upper and 
lower inlet quadrants, respectively. In figures 8 to 11, we present 

= A^pp for a large range of values of wave number k in the 

four quadrants where subscripts u and s refer to the unscarfed and scarfed 
inlets. A positive value of this quantity indicates that the conventional inlet 
radiates more acoustic power than the scarf inlet in that quadrant. From these 
figures, it is seen that for the range of frequencies considered, the scarf inlet 
radiates less acoustic power in the upper quadrant of the inlet and more 


20 log 
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power in the lower quadrant than the conventional inlet. This is clearly seen 
in Fig. 12 where we have plotted Figs. 8 and 11 on the same plot. It is seen 
that the radiated power in each of the inlet quadrants is a function of the wave 
number. Also there is an interesting symmetry with respect to the 0 dB line in 
this figure which indicates that the function of the scarf inlet is essentially to 
redirect the acoustic power from the upper quadrant to the lower quadrant. 
Figures 9 and 10 show the relative acoustic power radiation in the upper and 
the lower quadrants of the exhaust. There are both an increase and a decrease 
of acoustic power radiation of the duct with scarf inlet as compared to the 
conventional inlet depending on the wave number. There is no clear trend 
emerging from power radiated from the exhaust in our study. 

In Figures 8 to 1 1, we have shown the condition number of the matrix of the 
coefficients of a set of functions used in the boundary integral equation 
method. Large condition numbers indicate sensitivity to errors of the inver- 
sion of the matrix, except at isolated wave numbers, it appears that the condi- 
tion numbers are generally small and the precision of our calculations can be 
trusted. 

5.0 Concluding Remarks 

In this study, we have shown that the scarf inlet has the potential of redirect- 
ing the acoustic power into the quadrant with shorter inlet length. This behav- 
ior extends to a large range of wave numbers and appears to be also effective 
at forward speed. There are also changes in the exhaust radiation pattern 
which do not appear to follow a clear trend. These changes require further 
studies. We have also shown the usefulness of the boundary integral equation 
method code TBIEM2D used in our study. All calculations were performed 
on a PC at a short execution time. 


A Study of Noise Radiation from a Two Dimensional Scarf InletJuly 1, 1997 


S 




2-D Duct with Scarfed Inlet 
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Figure 8 




Figure 9 
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Figure 10 
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Figure 11 




Radiation from Inlet: A SPL = Unscarfed 



Figure 12 





